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PACS 03 . 65 . Ca - Formalism 

PACS 03.65.Ud - Entanglement and quantum nonlocality 
PACS 03 . 67 . Bg - Entanglement production and manipulation 
PACS 03 . 67 . Mn - Entanglement measures, witnesses etc. 

Abstract. - Although a precise description of microscopic physical problems requires a full quan- 
tum mechanical treatment, physical quantities are generally discussed in terms of classical vari- 
ables. One exception is quantum entanglement which apparently has no classical counterpart. We 
demonstrate here how quantum entanglement may be within the de Broglie-Bohm interpretation 
of quantum mechanics visualized in geometrical terms, giving new insight into this mysterious 
phenomenon and a language to describe it. On the basis of our analysis of the dynamics of a 
pair of qubits, quantum entanglement is linked to concurrent motion of angular momenta in the 
Bohmian space of hidden variables and to the average angle between these momenta. 
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Introduction. - The recent explosion in activity re- 
lated to entanglement is a consequence of a growing real- 
ization of its importance as a vital resource in quantum 
information through quantum teleportation jT], quantum 
cryptographic key distribution [2] and quantum computa- 
tion [3J[3] as well as its potential for enhanced quantum 
sensing through the engineering of highly entangled quan- 
tum states, beating the usual quantum limit [S]. Despite 
this utility, the concept of entanglement remains mysteri- 
ous and paradoxical within the usual quantum treatment 
and remains a focal point for discussing the foundations of 
quantum phenomena and related questions of non-locality 
[5], with renewed impetus from the recent developments 
in quantum information and quantum technology 

The connection between the classical and quantum ap- 
proaches is usually made through the correspondence prin- 
ciple although an alternative description based on the his- 
toric work of de Broglie [7] , and later developed by Bohm 
[8] gives equivalent results provided an effective non-local 
potential is introduced along with hidden variables. This 
approach has had a resurgence in interest recently re- 
emphasizing the virtues of the ontological description it 
provides [§]. 

With the development of the Bohm causal formulation 
of quantum mechanics [5], first for single spinless parti- 
cles and later for many-body fermionic or bosonic systems 
including fields [HUE], it was shown that the Bohmian 
approach is isomorphic with usual quantum mechanics re- 



garding observable predictions, but it additionally gives 
an interpretation of internal motion of particles in terms 
of classical coordinates, velocities and angular momenta. 
Although these (hidden) variables are not measurable, the 
approach is able to analyse individual processes in a way 
which goes beyond the standard Bohr interpretation of, for 
example, double slit experiments |12j . tunneling of parti- 
cles through barriers [13], Stern-Gerlach experiments [H] 
and the Aharonov-Bohm effect [15] . 

However, in the Bohmian approach quantum entangle- 
ment has not yet been analyzed quantitatively, although 
since Bell's construction of a hidden variable model for a 
single spin- 1/2 system (qubit) [1(3 [T7] the formal Bohm 
approach was consistently extended to a causal theory ca- 
pable of giving insight into Einstein-Podolsky-Rosen spin 
correlations in terms of well-defined individual particle 
trajectories with continuously variable spin vectors [TB] . 
Today spin- 1/2 systems can be treated causally in a non- 
relativistic formalism based on the Pauli equation |19j . 
by rigid rotor theory |20j . the hypersurfacc Bohm-Dirac 
model for entangled particles [21] or by the Clifford alge- 
bra approach to Schrodinger and Dirac particles [22] . 

Among the simplest quantum systems is a pair of qubits 
for which quantum entanglement can be quantified, for 
example, by the entanglement of formation Ep |23j . the 
asymptotic conversion rate to maximally entangled states 
from an ensemble of copies of a non-maximally entangled 
state [J]. The entanglement of formation can be related 
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to an associated quantity, concurrence C [21] ■ 
Here we consider a qubit pair in a pure state 

|*) = cos | |tl) + sin | ||t) , 



entangled 



(1) 



and for convenience we use spin-1/2 notation where |tl) 
corresponds to the first qubit in the "up" state, i.e., in the 
direction of the z-axis, and the second qubit in the "down" 
state. The generalization to more general forms of qubit 
pairs including mixed states is possible, but is not consid- 
ered here. Qubits are not restricted to real electron spins, 
but may be realized by any two state quantum system 
such as for example, entangled photon [25] , flux qubit in a 
superconducting ring [55], charge pseudo-spin of electron 
pairs in a double quantum dot [37]: nvm g qubits in quan- 
tum point contacts [28] or qubits in a composite system 

De Broglie-Bohm formalism for spin-1/2. - In 

order to demonstrate how quantum entanglement can be 
described ontologically we closely follow the approach in- 
troduced by Holland [TUH^U] where the starting point is 
the mapping between a quantum rigid rotor and a classi- 
cal spinning top in the presence of a quantum potential. 
The orientation of a rigid rotor is expressed by Euler an- 
gles £ = {a,/3,7} and, in the quantum approach, by the 
wave function "0(C)- 

Defining a differential operator M, whose components 
are the infinitesimal generators of the rotation group 
SO(3), the quantum Hamiltonian is given by H = M 2 /2J, 
where / is an auxiliary parameter (" moment of inertia" ) , 
with / — > in the final results. The wave function is ex- 
pressed as ip = lZe lS , where TZ(() and <S(£) are real func- 
tions. Bohmian space angular momentum is then given by 
a real three dimensional vector M = WIS. This relation 
is an analogue of a more familiar de Broglie's guidance 
equation for the velocity of a point-like particle with mass 
m treated in the Bohmian approach, mv = V5 [5]. 

The dynamics is determined from a Hamilton- Jacobi- 
type equation corresponding to the classical Hamiltonian 



H = 



M 2 

21 



Q = 



M 2 TZ 
2IK ' 



(2) 



where the quantum potential Q generates a quantum 
torque T = -iMQ, which rotates the angular momen- 
tum vector via the equation of motion dM./dt = T along 
the trajectory £(£). This is a counterpart of the Newton 
equation for the case of a free particle in the Bohm formu- 
lation given by mdv/dt = V(\7 2 lZ/2mlZ). The equation 
of the angular momentum motion simplifies to a set of first 
order non-linear differential equations where the solutions 
£(£) represent orbits in the configuration space, uniquely 
determined by the initial positions £(0), and the angular 
momentum emerges as M[£(£)]. 

Solutions for a single spin-1/2 rotor are known. Quan- 
tum basis eigenstates of M belong to the spin-1/2 sec- 
tor whereas the Bohmian description embodies the coex- 
istence of the motion of a rigid top, whose configuration 
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Figure 1: (a) An unentangled qubit pair - momenta are in- 
dependent and precess in opposite directions around the z- 
axis. (b) Typical motion of momenta in a partially entangled 
state, i? = 7r/5 and tp = 0, with initial values for the rotors 
C(0) = {ai,£i,7i,a2,/3 2 ,72} = {2,5,0, |, 0,0}. (c) Spire- 
graphic motion of M(t) = Mi + M2 projected onto the xy- 
plane (parameters as in (b). (d) Trajectories of Mi, 2 shown 
in the xy-coordinate system rotating with M(£). In this coor- 
dinate frame the trajectories are closed and periodic with the 
period r(C). 



space is SO(3) and a guiding wave, whose spin configura- 
tion space is SU(2) [3D]. The pilot wave is given in terms 
of the Wigner D-matrices. For the "spin up" state, for 
example, 7£f(C) °c cos ^ and Sf (£) = —(/? + r y)/2. 

From the equations of the motion follows a time depen- 
dent angular momentum M: "spin up" and "spin down" 
means that M, starting from some initial direction which 
is a function of the initial value C(0), processes around the 
2-axis counterclockwise and clockwise, respectively, with 
constant z-axis projection M z = ±i. The angular mo- 
mentum magnitude is constant and can reach any value 
|M| > |, depending on the initial choice of £(£). Av- 
eraging over the initial values yields a quantum equilib- 
rium ensemble averaged angular momentum (M)b, which 
is time independent and identical to the quantum mechan- 
ical expectation value of the spin operator (S) and also to 
the time average of Bohm angular momentum (M)t for 
each particular ensemble representative £. In Fig. 1(a) is 
shown an example of two unentangled qubits in the state 
|ti). The solution is an independent motion of "spin up" 
and "spin down" vectors Mi and M2. 

For a two qubit state is the guiding function ip(Q ex- 
pressed by six variables C = {CijC2}- It should be noted 
that even for non-interacting, but entangled qubits the 
corresponding equations of the motion are coupled by the 
quantum potential Q(() which incorporates their interac- 
tion in the corresponding Bohmian two-particle Hamilto- 
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nian. The solutions for each of the time-dependent angular 
momentum vectors Mi 2 are determined by six common 
initial values £(0). 

The probability distributions and the quantum 
entanglement. - Let us now discuss results for a qubit 
pair in the state parametrized by equation ((T|). The total 
angular momentum projection Miz+M^z is zero while the 
angular momenta due to the action of the non-local quan- 
tum potential (i.e., interaction) and the corresponding 
quantum torques exhibit a complex precessional motion, 
as illustrated for a particular choice of £(0) in Fig. 1(b). In 
general, only a few special cases have been explored so far. 
For example, for a fully entangled singlet state, -d = it/2 
and tp = ir, the total momentum M = Mi + M2 vanishes 
for each £ in accordance with the usual imagery [TU] . 

In this letter we concentrate on dynamical properties 
of the angular momenta relevant to the analysis of quan- 
tum entanglement while various other spin-spin correla- 
tion functions and the corresponding probability distribu- 
tions will be presented elsewhere [3D]. We computed tra- 
jectories Mi,2[C(*)] covering the full configuration space 
with ~ 10 6 initial values £(0) per \ty), i.e., for a particu- 
lar choice of d and tp. Although these trajectories exhibit 
extremely rich variety, some common properties can be 
outlined, (i) The quasi-periodic motion appears chaotic 
and, except in special cases, the projections of the total 
momentum M onto the a;y-plane winds around the origin 
an infinite number of times in a spirographic manner |31] , 
forming a dense annulus limited by fixed outer and inner 
radii [Fig. 1(c)]. (ii) The curve corresponding to relative 
momentum M2 — Mi is closed and periodic if plotted in 
the reference frame rotating synchronously with M around 
the z-axis [Fig. 1(d)]. 

These results clearly show that the entanglement prop- 
erties of such a qubit pair are reflected in the dynam- 
ics of the azimuthal angles 0i [C(*)] an d 4>2[((t)] of an- 
gular momenta as follows. First, the ensemble average 
difference of azimuthal angles (j>[C(t)} — <f>2 — 4>i IS time 
independent and given by (4>)b = <P, where the aver- 
age is defined by (f) B = J/(C)^ 2 (C)dC with dC = 
n?=i sin ctidcxidPid'yi. The corresponding probability dis- 
tribution dP{<j>)/d<t> = fS[<f>- 0(C)]7e 2 (C)dC [22 is con- 
stant for unentanglcd qubits and becomes progressively 
peaked at tp for increasing entanglement, as presented in 
Fig. 2, culminating in precession of angular momenta at 
equal relative angle 4>[((t)} = tp for all f consistent with 
perfect entanglement. The shape of the distribution is 
independent of tp. 

The next quantity of interest is the probability distri- 
bution t£P(cos((/>— tp)) / d cos((f>— ip), also presented in Fig. 2 
and, in particularly the average cosine which appears to be 
an excellent measure of entanglement, ranging from zero 
for an unentangled state to unity for a maximally entan- 
gled state. One can readily prove the exact expressions 
(cos</>)b = Cb cos ip and (sin^B = Casing or, equiva- 




cos((|)-cp) 

Figure 2: Probability distribution for relative angle <f> of the in- 
plane projections of momenta (see Fig. 1(a)), dP{cj>)/d<j> (blue, 
upper scale) shown relative to the average angle {(f) b = f 
for various degrees of entanglement. In the low entanglement 
regime the distribution is flat, becoming peaked at tp in the 
fully entangled state. The corresponding probability distribu- 
tions dP(cos((j>—(p))/dcos((j)—ip) (red, lower scale). The aver- 
age of (cos(0— tp)) b = Cb{"&) is marked with thick vertical lines 
and shown also in the inset with bullets together with concur- 
rence sini? (dashed line). Small oscillations are due to finite 
size effects. 
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where Cb is dependent only of d and A cos </>, A sin </> are 
standard deviations from average cosine and sine, respec- 
tively. Note that (sin(</> — <p))b = 0. Similar formulae can 
be derived also by appropriately defined cosine and sine 
operators in standard quantum mechanics formalism |33j . 
A higher degree of entanglement can thus be visualized as 
a highly correlated distribution of angular momenta mak- 
ing azimuthal angles difference close to tp, with suppressed 
fluctuations for progressively increasing entanglement. 

In the Bohmian picture of entangled qubit pairs the 
quantity Cb{$) (inset to Fig. 2) plays the role of Woot- 
ters concurrence given by C{d) = 2|(4 , |S' 1 I ~S' 2 ~| , I')| = (sin i9] 
|24j , where 5*^2 are spin-ladder operators for qubits 1 and 
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2, respectively [3~4"Il3"5] . High Cb(i9) —> 1 signals weak an- 
gle fluctuations whereas for Cb{&) — > 0, entanglement is 
suppressed and progresively destroyed when the standard 
deviation of the cosine (or sine) is comparable to its aver- 
age. 

Concurrence is related to quantum mechanical expec- 
tation values and Cb to ensemble averages. A natural 
question arises: Is there some imprint of entanglement, 
not only in the full ensemble, but also in each particular 
representative £? Standard quantum mechanics does not 
discuss such questions, while in the Bohmian approach one 
can classify the ensemble further. 

Each representative £ will in general lead to a different 
time average of cos (/>([( (t)]— tp), defined by 



)r .I 



cos(<f>[((t)}-<p)dt, 



(7) 



where t(£) is the period corresponding to the trajectory 
((t). Time averaged cos(</>[£ (t)]—ip) is distributed accord- 
ing to the probability distribution dP((cos ^)r)/^(cos </>}t- 
In Fig. 3 is presented the distribution for several represen- 
tative values of The distributions in the strong entan- 
glement regime, i? ~ n/2, are peaked at (cos</>)t ~ 1, 
similar to results in Fig. 2 but strikingly different in the 
weak entanglement regime. This is understandable since 
for weakly entangled qubit pairs the angular momenta pro- 
cess almost independently with a vanishing time average 
cos(0 — ip) for every £, leading to the zero-peak in the dis- 
tribution. This is contrary to the fixed time distribution 
which is flat because of the randomly distributed angles. 

It should be noted that the ensemble representing two 
qubits is non-ergodic, i.e., ensemble averages (/)b do not 
generally equal time averages (/}t- For the present case of 
nonintcracting (but entangled) qubits {f)s is independent 
of time (as it should be) and {{/)t)b = where the 

Bohmian ensemble average of (/)t is evaluated using the 
corresponding initial values £(0). 

The probability distribution for the cosine time aver- 
age exhibits a distinctive tripartite structure. We found 
that the discriminating property of ensemble representa- 
tives is the relative direction of angular momenta preces- 
sion. In the low entanglement regime the xy-plane pro- 
jection of momenta Mi and M2 precess mainly in oppo- 
site directions. An extreme case is an unentangled state, 
Fig. 1(a). In general, momentum pairs move part time in 
the same and part time in the opposite direction. We clas- 
sify representatives that always precess in the same direc- 
tion as "concurrent" movers whereas those which always 
precess in the opposite direction are classified as "anti- 
concurrent" . 

To be specific, we introduce a measure Cq to distin- 
guish different trajectories ((t) according to their "con- 
currency" , 
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Figure 3: Time averaged cos(<j>— ip) probability distributions 
for various degrees of qubit pair entanglement. The concurrent 
(Cf = 1) and anti-concurrent (Cq = — 1) fractions are indicated 
by red and blue shading, respectively. Thick vertical lines rep- 
resent distribution averages (values are identical to those in 
Fig. 2, i.e., C B ). 



At each moment the angular momenta for a given trajec- 
tory ((t) precess either in the same or in the opposite direc- 
tion. One can thus visualize the concurrency as a measure 
of the share of the time that both angular momenta move 
in the same direction. For example, = ±1 for perfectly 
concurrent and anti-concurrent movers, respectively, and 
Cq > for trajectories where angular momenta move con- 
currently more than half of the time for some members of 
the ensemble. 

This is shown in Fig. 3 where distributions for con- 
current and anti-concurrent movers are shaded red and 
blue, respectively. In the low-entanglement regime anti- 
concurrent movers dominate whereas the distribution of 
concurrent movers progressively dominates as entangle- 
ment increases. Remarkably, only a minority of repre- 
sentatives move both concurrently and anti-concurrently, 
part of the time (\C^\ < 1). 

Concurrency should not be confused with concurrence 
which is given by the ensemble average of cos <f> while, on 
the other hand, concurrency relates to spin-precession of 
a particular pair in the Bohmian ensemble and has no 
direct quantum analogue. The probability distribution of 
concurrency consists of a discrete and a continuum part, 



dP(C c ) 



P + 6(C ( - I) + P-5(Cs + l) + p(C c ), (9) 



dt 



dl 



(8) 



where P± is the probability that the concurrency is ex- 
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Figure 4: (a) The histogram d 2 P/d(cos </>)tcK< extracted from 
~ 10 6 trajectories for ■& — ir/5 [32]. The samples are identical 
to those in Fig. 2, while a larger bin size is applied here 
in order to enable the presentation of a discrete/continuous 
type of the distribution in a single plot. Thick vertical line 
represents the average point {Cb, (Q)s}- (b) Continuous part 
p(C() vs. concurrency for various The distribution Eq. ([9j) 
is normalized, thus J p(C^)dC<; = 1 — (P+ + P-). 



actly ±1, respectively, and p(Cf) is a continuous function 
[36], shown for various # in Fig. 4(b), for which motion is 
sometimes concurrent and sometimes anti-concurrent as t 
changes. Hence the probability density p{C^) tends to zero 
as |Cf | — > 1, which strictly separates (anti-)concurrent mo- 
tions from the partially-concurrent regime. For a typical 
case of a partially entangled qubit pair with § = tt/5 is 
in Fig. 4(a) shown the histogram probability distribution 
d 2 P/d(cos<f))TdC^. Coarse binning is applied in order to 
emphasize the middle regime \C^\ < 1 together with (anti-) 
concurrent parts. 

Let us emphasize, that concurrent motion is inherent 
to an individual representative and their share P+ in the 
ensemble can be considered as a suitable measure of en- 
tanglement for Bohmian ensembles. Such a clear visual- 
ization of entanglement cannot be deduced from quantum 
mechanics where only averages of operators are accessi- 
ble. Quantifying entanglement is important in applica- 
tions which require maximally entangled (Bell) pairs to 
be shared by two remote parties. In reality qubit pairs 
in general are not perfect Bell states and entanglement 
distillation is required in order to extract such pairs [25] . 
The entanglement of formation Ep is an upper bound to 
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Figure 5: (a) P+ line represents the probability that momenta 
move concurrently (shaded red), and above 1 — P_ is the proba- 
bility that motion is anti-concurrent (blue shaded region) . The 
Pa line denotes a share of trajectories which move concurrently 
more than half of the time (Po = P+ + J Q p{Cq )dCf ). Entangle- 
ment of formation Ef is shown with a dashed line, (b) Entan- 
glement of formation vs. the probability for configuration 
in \9) (dashed). Bullets represent the probability for concur- 
rent motion - the Bohmian counterpart of Ef- The full line 
represents the distillation yield using Procrustean method. The 
dotted line is the yield using the Schmidt projection method 
with 22 input pairs. 



the average number of Bell pairs that can be extracted or 
distilled from a set of copies of an entangled state, using 
only local operations and classical communication. In the 
Bohmian approach the probability for concurrent motion 
P+ may be considered a counterpart of Ep. This prob- 
ability is shown in Fig. 5(a) together with Ep (dashed 
line). The probability for anti-concurrent motion P_ cor- 
responds to the region above the 1 — P_ line. 

In Fig. 5(b) the distillation yield of Bell pairs obtained 
by the Procrustean method [23] is plotted as a func- 
tion of cos 2 $/2 for one qubit pair (straight full line) 
together with the upper limit, the entanglement of for- 
mation (dashed line). The probability for concurrent 
movers in the Bohmian picture (red line) is above both 
the Procrustean method and the yield of Schmidt projec- 
tion methods when applied with n = 22 input pairs |23| . 
This result suggests that in the Bohmian picture "pre- 
formed entangled" pairs exist and their extraction would 
represent a quantum distillation protocol with high yields 
obtained from a single qubit pair. 

Conclusion. - For the sake of completeness let us dis- 
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cuss also the class of states 



l*)=cos^|tt) + ^sin^|U) 



(10) 



The formalism can be applied in a similar manner as for 
the case of \^>) while the main distinction is that for states 
Eq. (fTU)) the sum of azimuthal angles, <f> = <f>\ + <p2 . is im- 
portant, not the difference <f>. Although the trajectories 
with equal initial values, but different quantum states |\I>) 
d'Jj}), are altogether different, the probability distribu- 
tions take the same form after the appropriate substitu- 
tion (f> <-> 4>. For example, (4>)b = f and, furthermore, the 
sign of concurrency is reversed, i.e., Q = — Cf. 

In brief summary the main findings are as follows. In the 
usual quantum mechanical approach, two main measures 
are used to quantify entanglement, namely concurrence 
and the entanglement of formation. These rather abstract 
concepts are related but rather difficult, if not impossi- 
ble, to visualise. Conversely, in the Bohmian approach it 
is possible to define analogous measures that are quanti- 
tatively similar (though not identical) to concurrence and 
entanglement of formation but are quite distinct from each 
other and may be directly related to intuitively appeal- 
ing geometric interpretations. In the Bohmian interpreta- 
tion, the angular momentum vectors of the two particles 
(qubits) process in a well-defined way with some initial 
probability distribution. The Bohmian concurrence is sim- 
ply the average mutual cosine of the xy-plane projection 
of the two angular momentum vectors and is a property 
of the whole ensemble. In standard quantum mechanics 
this corresponds to an expectation value. On the other 
hand, the Bohmian entanglement of formation depends 
on a property - concurrency - of each representative of 
the ensemble which either has "concurrent motion" (the 
angular momentum vectors precess in unison) or they do 
not. The share of representatives with this motion is the 
entanglement of formation. Apart from this simple and 
appealing underlying picture, it may also have practical 
relevance since it suggests that distillation protocols could 
have a very high yield already for one single qubit pair. 
This ontological result is entirely missed (has no meaning) 
in the usual quantum mechanical approach. Our results 
also give rise to a challenging question: Can some quan- 
tity identical, or analogous, to concurrency also manifest 
itself in other spin- 1/2 formalisms? 
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